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I. INTRODUCTION 

Nonlinear magneto-optical rotation (NMOR) , or light- 
power-dependent rotation of optical polarization due to 
resonant interaction with an atomic medium in the pres- 
ence of a magnetic field B, has applications ranging from 
fundamental symmetry tests to magnetometry 1] . With 
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FIG. 1: Experimental (dots) and dispersive-Lorentzian least- 
squares fit (line) magnetic-field dependence of NMOR in a 
10-cm-diameter paraffin-coated 85 Rb-vapor cell, obtained as 
described in Ref. ||. Linearly polarized laser light is tuned to 
the high-frequency side of the D2 F = 3 —* F' transition, at 
which maximum rotation occurs. The light intensity is ~50 
/jWcm" 2 , and the beam diameter is ~2 mm. The tempera- 
ture of the cell is ~19°C corresponding to a vapor density of 
~4 x 10 9 cm- 3 . 
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NMOR due to the evolution of ground-state atomic po- 
larization , optical rotation is proportional to the mag- 
netic field for small fields, but falls off when the Larmor 
frequency = g^B (g is the gyromagnetic ratio, /i 
is the Bohr magneton, and we set K = 1 throughout) 
becomes larger than half of the atomic polarization re- 
laxation rate 7 (Fig. ^) . Atomic polarization relaxation 
rates as low as 7 ~ 2ir x 1 Hz can be achieved for alkali 
atoms contained in paraffin-coated vapor cells Q| , corre- 
sponding to magnetic field widths of approximately 1 [iG 
[» and high magnetometric sensitivity (~3 pGHz -1 / 2 
|3) to small fields. 

With a traditional NMOR magnetometer, the high 
small-field sensitivity comes at the expense of a limited 
dynamic range. Since many applications (such as mea- 
surement of geomagnetic fields or magnetic fields in space 
0) require high sensitivity at magnetic fields on the or- 
der of a Gauss, a method to extend the magnetometer's 
dynamic range is needed, ft was recently demonstrated 
0, that when frequency- modulated light is used to in- 
duce and detect nonlinear magneto-optical rotation (FM 
NMOR), the narrow features in the magnetic- field de- 
pendence of optical rotation normally centered at B = 
can be translated to much larger magnetic fields. In this 
setup (Fig. [2Jl, the light frequency is modulated at fre- 
quency firm an d the time-dependent optical rotation is 
measured at a harmonic of this frequency. Narrow fea- 
tures appear, centered at Larmor frequencies that are 
integer multiples of f2 m , allowing the dynamic range of 
the magnetometer to extend well beyond the Earth field. 

Light-frequency modulation has been previously ap- 
plied to measurements of linear magneto-optical rota- 
tion and parity-violating optical rotation [SJ, 1 1 0I [ in or- 
der to produce a time-dependent optical rotation signal 
without introducing additional optical elements (such as 
a Faraday modulator) between the polarizer and ana- 
lyzer. Optical pumping with frequency-modulated light 
has been applied to magnetometry with 4 He E3 
and Cs [14[; in these experiments transmission, rather 
than optical rotation, was monitored. In the latter work 
with Cs, the modulation index (the ratio of modulation 
depth to modulation frequency) is on the order of unity, 
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FIG. 2: Simplified schematic of the apparatus used to detect 
FM NMOR signals. A paraffin-coated cell containing Rb va- 
por is placed inside a balanced polarimeter (a polarizer and an 
analyzer oriented at ~45° with respect to each other). The 
frequency of the laser is modulated with an amplitude of a 
few dozen MHz. The lock-in amplifier is used to detect the 
components of optical rotation oscillating both in phase and 
7r/2 out of phase with the frequency modulation. 



in contrast to the much larger index in the work de- 
scribed here, allowing interpretation of the process in 
terms of the A- or coherent-population-trapping reso- 
nances. This regime has also been explored in Rb using 
modulation of the magnetic field, rather than the light 
field [l5(. The closely related method of modulation of 
light intensity (synchronous optical pumping) predates 
the frequency-modulation technique |lfi| . Also employing 
light-intensity modulation is the so-called quantum beat 
resonance technique fl?! used, for example, for measur- 
ing the Lande factors of molecular ground states (see Ref. 
[lj| and references therein). Intensity modulation was 
recently used in experiments that put an upper limit on 
the (parity- and time-reversal-violating) electric dipolc 
moment of 199 Hg (Refs. 0,113] and references therein). 

A quantitative theory of FM NMOR would be of use 
in the study and application of the technique. As a first 
step towards a complete theory, we present here a per- 
turbative calculation for a J = 1 -> J' = atomic tran- 
sition that takes into account Doppler broadening and 
averaging due to velocity-changing collisions. We begin 
the discussion in Sec. In] by comparing experimental FM 
NMOR magnetic-field-dependence data obtained with a 
paraffin-coated 87 Rb-vapor cell to the predictions of the 
calculation (described in Sec. lIII|) . We find that the sim- 
plified model still reproduces the salient features of the 
observed signals, indicating that the magnetic-field de- 
pendence of FM NMOR at low light power is not strongly 
dependent on power or angular momentum. As discussed 
in Sec. IIVI the description of the saturation behavior and 
spectrum of FM NMOR in a system like Rb, on the other 
hand, will require a more complete theory. 
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FIG. 3: Measured (left column) and calculated (right col- 
umn) in-phase (top row) and quadrature (bottom row) first- 
harmonic amplitudes of FM NMOR. The experimental sig- 
nals, plotted as a function of magnetic field B applied along 
the light propagation direction, are obtained with light tuned 
to the wing of the F = 2 — > F' = 1 absorption line of the 
87 Rb Dl spectrum. The laser power is 15 /iW, beam diame- 
ter is ~2 mm, Q m — 2ir x 1 kHz, and modulation amplitude 
is 27r x 220 MHz. All resonances have widths (~1 /^G) corre- 
sponding to the rate of atomic polarization relaxation in the 
paraffin-coated cell. The normalized calculated signals [Eq. 
1581 . for aj = l— > J' — transition, are plotted as a func- 
tion of normalized Larmor frequency S7i/Sl m . For these plots, 
the parameters Aq/Fd = 0.7, fi m /7 = 500, and Aj/Ao = 1 
(described in Sec. II1II are chosen to match the experimental 
parameters given above. 



II. EXPERIMENTAL DATA AND 
COMPARISON WITH THEORY 

Figures |3] and 0] show first- and second-harmonic data, 
respectively, obtained from an FM NMOR magnetome- 
ter with light tuned near the Dl line of rubidium in the 
manner described in Ref. 0, , along with the predicted 
signals for a J = 1 ^ J' = transition obtained from 
the theory described in Sec. IIIII with parameters match- 
ing those of the experimental data. The calculation for 
the simpler system reproduces many of the qualitative 
aspects of the experimental data for Rb. The features at 
the center of the in-phase plots of Figs. |3 and 0] are the 
zero-field resonances, analogous to the one shown in Fig. 
n (The background linear slope seen in the in-phase 
signals is also a zero-field resonance, due to the "tran- 
sit effect" It is modelled in the theory by an extra 
term analogous to the others with the isotropic relax- 
ation rate 7 equal to the transit rate of atoms through the 
laser beam.) In addition to these features, there appear 
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FIG. 4: Measured and calculated second-harmonic ampli- 
tudes of FM NMOR. See caption to Fig. For the experi- 
mental signals, light is tuned to the center of the F — 2 — > 
F' — 1 absorption line of the 87 Rb Dl spectrum, the laser 
power is 15 fiW, beam diameter is ~2.5 mm, Q m = 27rx 1 kHz, 
and modulation amplitude is 2-7T x 440 MHz. The parameters 
for the theoretical signals are Aq/Td = 1.4, Q. m /"i = 500, and 
Ai/Ao = 0.2. 
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FIG. 5: A J s = 1 -> J e = atomic transition of frequency 
ujo- The lower sublevels are split by the Larmor frequency Ql- 
The arrows indicate the interaction with light of frequency ui 
polarized perpendicular to the quantization axis. The upper 
state spontaneously decays a rate 70. 



There are additional features, centered at |fii/ri m | = 
1/4, just barely visible in the experimental plots of Figs. 
13 These features, which become more prominent at 
higher light power [8j, are due to the optical pumping, 
precession, and detection of the hexadecapole moment. 
These resonances are not described by the current the- 
ory, because the presence of the hexadecapole moment 
requires ground-state angular momentum J > 2 and 
second-order light interactions. A quantitative descrip- 
tion of these resonances is among the goals for an ex- 
panded theory. 



III. THEORY 



new features centered at magnetic field values at which 
|f2[,/f2 OT | = 1/2 and 1. For the first-harmonic signal, the 
former are larger, whereas for the second-harmonic, the 
latter are; this is primarily a result of the different light 
detunings used in the two measurements. For these new 
resonances, there are both dispersively shaped in-phase 
signals and 7r/2 out of phase (quadrature) components 
peaked at the centers of these resonances. The resonances 
occur when the optical pumping rate, which is periodic 
with frequency fl m due to the laser frequency modula- 
tion, is synchronized with Larmor precession, which for 
an aligned state has periodicity at frequency 2Ql as a 
result of the state's rotational symmetry. This results 
in the atomic medium being optically pumped into an 
aligned rotating state, modulating the optical properties 
of the medium at 2Ql. The aligned atoms produce max- 
imum optical rotation when the alignment axis is at 7r/4 
to the direction of the light polarization and no rota- 
tion when the axis is along the light polarization. Thus, 
on resonance, there is no in-phase signal and maximum 
quadrature signal. The relative sizes and signs of the 
features in the magnetic-field dependence, largely deter- 
mined by the ratio of the modulation width Ao to the 
Dopplcr width Yd fScc. are well reproduced by the 
theory. The theory also exhibits the expected linear light- 
power dependence of the optical rotation amplitude as 
observed in experiments at low power 0, 0] ■ 



A. Introduction 

The goals for a complete theory of FM NMOR are 
outlined in Sec. lIVI As a first step towards such a theory, 
we calculate here the optical rotation due to interaction 
of frequency-modulated light with a J g = 1 — > J e = 
atomic transition (Fig. [5}, where the subscripts g and e 
indicate the ground and excited states, respectively. We 
will assume that the light power is low enough that no 
optical pumping saturation occurs. 

We begin by calculating the time-dependent atomic 
ground-state coherence of a Doppler-free system. Us- 
ing the magnetic-field-atom and light-atom interaction 
Hamiltonians (under the rotating wave approximation) 
we write the density-matrix evolution equations. Under 
the low-light-power approximation, an expression for the 
ground-state atomic coherence can be written as a time 
integral. We convert the integral to a sum over harmonics 
of the modulation frequency by expanding the integrand 
as a series. This form is convenient for this calculation 
because the optical rotation signal is measured by lock- 
in detection. The expression for the Doppler-free case 
is then modified to take into account Doppler broaden- 
ing and velocity averaging due to collisions. Atoms in 
an antirelaxation-coated vapor cell collide with the cell 
walls in between interactions with the light beam, pre- 
serving their polarization but randomizing their veloci- 
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ties. In the low-light-power case, we can account for this 
by first calculating the effect of optical pumping assum- 
ing no collisions, and then averaging the density matrix 
over atomic velocity. Note that in this case we assume 
that optical pumping is unsaturated not only for the res- 
onant velocity group, but also when atomic polarization 
is averaged over the velocity distribution and cell volume. 

Using the wave equation, we find an expression for the 
time-dependent optical rotation in terms of the atomic 
ground-state coherence of a given atomic velocity group. 
This rotation is then integrated over time and atomic 
velocity to obtain an expression for the signal at a given 
harmonic measured by the lock-in detector. 



B. The Hamiltonian 

The total Hamiltonian H is the sum of the unperturbed 
Hamiltonian Hq, the light -atom-interaction Hamiltonian 
Hi, and the magnetic-field-atom-interaction Hamilto- 
nian Hb- Using the basis states |£JM), where £ rep- 
resents additional quantum numbers, denoted by 



|£eJe,0) = 



(1) 



the unperturbed Hamiltonian Hq is given by 



Hn = 



/0 



Vo o o 



(2) 



where loq is the transition frequency (again, we set h = 1 
throughout). 

An x-polarized optical electric field E is written as 



E = Eq cos (ut) e x . 



(3) 



where Eq is the electric field amplitude and to is the fre- 
quency (modulated as u> = lji — Ao cos fl rn t, where u>i 
is the laser carrier frequency and Ao is the modulation 
amplitude). We assume that the atomic medium is op- 
tically thin, so that we can neglect the change in light 
polarization and intensity inside the medium when cal- 
culating the state of the medium. The light-atom inter- 
action Hamiltonian is given by 



Hi = -E • d 

= —Eq cos (cut) d x 

= — \=Eq cos (orf) (d-i 
v2 



(4) 



d+i) , 



where d is the dipole operator. According to the Wigner- 
Eckart theorem, components of a tensor operator T Kq are 
related to the reduced matrix element (S,J\\T K \\^J') by 

E2 



^JM\T m \i'J'M') = {-l) J - M m\T K U'J') 



J K J' 

—M q M' 
(5) 

Thus the matrix elements of d+i and g?_i for this transi- 
tion can be written 



{ijM\d ±1 \i'j'M') = 

f (^JeWdUgJg) for £ = £ e , J = J e ,M = 



,/ I J' 
—M ±1 M' 



. and £' = £,, J' = J g ,M' = T l, 

X \ (CgJgWdUeJe) for f = J = J g , M = ±1 

x ° and £' = J' = J e ,M' = 0, 







in all other cases. 



(6) 



Reduced matrix elements with different ordering of states 
are related by [2l| 



mT K \\ej') = (-i) J - r (CJ'mujy , 

and since the reduced dipole matrix element is real, 

(Z g J g \\d\\ZeJe) = -(ZeJe\\d\\t B J g ). 

Thus Hi is given in matrix form by 



Hi = 2Q cos tot 



-IN 



1 

\—l I 



(7) 



(8) 



(9) 



where f2 = (£ g J g ||d||^ e J e ) Eq/(2^/6) is (apart from a nu- 
merical factor of order unity) the optical Rabi frequency. 

The magnetic field interaction Hamiltonian Hb for a 
z-directed magnetic field B is given by 



H B = -H-B 
= gfio J • B 
= gfi J z B 

-1 0\ 

10 



(10) 



0/ 



where J7l is the Larmor frequency as defined in Sec. [I] 
Thus, the total Hamiltonian is given by 

H = Hq + Hi + Hb 

-Cl L -2ficoswA 

ooo I ( n ) 

n L 2Qcoscut 
-2Qcosujt 2Qcosujt ujq 
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C. Rotating- wave approximation 

We now use the rotating-wave approximation in order 
to remove the optical-frequency time dependence from 
the Hamiltonian. We first transform into the frame ro- 
tating at the optical frequency by means of the unitary 
transformation operator U(t) = exp(—iH't), where 



H' 



'0 ON 





v0 ujj 



(12) 



is the unperturbed Hamiltonian Hq with loq replaced by 
uj. It is straightforward to show that under this transfor- 
mation the Hamiltonian in the rotating frame is given by 



H = U-\t) H(t) U(t) - i U-\t) — U(t) 



o 
o 



o 



n L 



o nh+e 2 



-fi (l + e 


fi(l + e- 

UJq — LU 



I 

(13) 



2iut 



where we have used costot = (e~ lult + e lwt ) / 2. Averag- 
ing over an optical cycle to remove far-off-resonant terms 
(the rotating wave approximation), we have 



H 



-n L o o 



OSli 

-q o n 



(14) 



where A = lu — loq is the (time-dependent) optical detun- 
ing. 



D. Relaxation and repopulation 

We assume that the upper state spontaneously decays 
with a rate 70, and that the ground state relaxes with a 
rate 7, due to the exit of atoms from the light beam, in 
the case of the "transit" effect, or collisions with other 
atoms or the cell wall in the case of the "wall-induced 
Ramsey effect" pj. (Additional upper-state relaxation 
processes can be neglected in comparison with the spon- 
taneous decay rate.) This relaxation is described by the 
matrix T, given by 



/7 0' 

7 

7 

\0 7o , 



(15) 



The simplest model of ground state relaxation is used. 
The effects of collisional dephasing could be included by 
adding off-diagonal terms if a more realistic model is de- 
sired. In order to conserve the number of atoms, the 



ground state must be replenished at the same rate at 
which it relaxes. This is described by the repopulation 
matrix A: 
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/7 0\ 

7 

7 

Vo 0/ 



(16) 



where N is the atomic density. We ignore repopulation 
due to spontaneous decay since the calculation is per- 
formed in the low-light-power limit (fl 2 <C 707)- 



E. Density-matrix evolution equations 

The evolution of the density matrix p (defined so that 
Tr p — N) is given by the Liouville equation 



1 



{r,p} + A, 



(17) 



where the square brackets denote the commutator and 
the curly brackets the anticommutator. The M = 
ground-state sublevel does not couple to the light, and 
can be ignored. Using |1) and |2) to denote the ground- 
state M = —1 and +1 sublevels, respectively, and |3) to 
denote the upper state, and assuming that 7 <C 70, the 
evolution equations for the atomic coherences obtained 
from Eq. l(T7|) are 

P31 = - [7o/2 + i (Ql - A)] p 3 i + ifl (pn - p 33 - P21) , 

(18a) 

P23 = - [7o/2 + i (^L + A)] p 23 + iVL (p 22 - P33 ~ P2l) , 

(18b) 

P21 = - (7 + 2ifl L ) P21 - i£l (p3i + P23) ■ (18c) 

We can assume that in the low-light-power limit the pop- 
ulations pn,22,33 are essentially unperturbed by the light 
(pn,22 — N/3, pz3 -C N). We can also assume that, 
neglecting transient terms, the optical coherences P3i,23 
are slowly varying (any time dependence would be due to 
modulation of the light frequency, which will always be 
done at a rate much less than 70; thus /93i,23 "C 70/031,23)- 
Using these assumptions, the evolution equations for the 
atomic coherences [Eqs. (|18aH18cjl ] become 

~ - [70/2 + i (Q L - A)] pa + in (N/3 - p 2 i) , (19a) 
~ - [70/2 + % (fi £ + A)] P23 + ifi (N/3 - P 2i), (19b) 
P21 ~ - (7 + 2ifl L ) P21 - in (p 3 i + P23) ■ (19c) 

These equations can be used to solve for the optical and 
ground-state coherences. 



Calculation of the optical and ground-state 
coherences 



The expression for optical rotation fSec. lIHHjl is writ- 
ten in terms of the optical coherences P3i,23- We will now 
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relate the optical coherences to the ground-state coher- 
ence P21 and find an expression for p 2 i as a sum over har- 
monics of the light detuning modulation frequency Q m . 
This form is convenient because the signal is measured 
at harmonics of this frequency. 

Solving Eqs. (|19a|l and ill9bl) for p 3 i and P23 in terms 
of P21 , we obtain 



P31 

P23 



n(N/3-p 21 ) 
Q L -A- i 7o /2 

n (N/3 - p 21 ) 
n L + A- ijo/2 



(20) 



In order to solve for p 2 i, we make the substitution p 2 i 



t%\ e 



-(2tn L +7)i 



Eq. Pc|l : 



f 2 i^~in(p 31 +p 23 )e^ n ^~>\ (21) 
or, integrating (assuming that r 2 \ = p 2 \ = at t = 0), 



rai ~ -ifi / (P31 + P23) e (2 ^+^dT, (22) 
Jo 



so, substituting back, 



in / (A3i+P2 3 )e- (24n ^ )(t - T W. (23) 



P21 



The expressions for the optical coherences [Eqs. <|20|) ] are 
then substituted into the expression for the ground-state 
coherence [Eq. (|23p]. Assuming that the light power is 
low (f2 <C 70) allows us to neglect second-order terms. 
We also assume that the level shift induced by the mag- 
netic field is smaller than the natural line width, i.e. 

<C 70 • (For the D- lines of rubidium used in the ex- 
periment, this assumption holds for magnetic fields up 
to the earth-field range.) The ground-state coherence is 
then given by 



1 



1 



Q 2 N 



Q L - A - ijo/2 Q L + A- i la /2 

X e -(2^+7)(*-r) dr 

'70/2 -iflL 2iQ L A 2 \ 



V A2 +^o74 (A 2 + 7o 2 /4)7 

X e -( 2 ^+7)(*-) dr 



~ n 2 N [( 7o /2 - m L ) h(t) + 2in L i 2 (t)} 



where the integral 1% has been defined by 



(24) 



-(2in £ + 7 )(t-T) rfr 



A 2 (L> -cosfi m r) 2 + 7 2 /4 
= K/4) _1 f m m r) e-^+^Ur, 



(25) 



and I 2 has been defined by 



h{t) = 



Aq (I?o — cos f2 m T) c?t 



A 2 (C -cosri m r) 2 +7 2 /4 



= -A 



2 9h (t) 
cA 2 ' 



-(2in £ +7)(t-T) 



(26) 



Here we have substituted for A the expression for the 
light- frequency modulation A = Ao (Dq — cos f2 m r), 
where the dimensionless average detuning parameter Do 
is defined by Dq — A;/Ao, where A/ = loi — luq. The 
lineshape factor f(x) is defined by 



fix) 



7oV4 



A 2 (D -cos2;) 2 +7o/4' 



Expanding this function as a series of harmonics, 

00 

f{x)= *ne mx , 

71— — OO 

the coefficients a n are given by 

1 r 

a„ = — / fix) cosnxdx. 
2W_/ 



(27) 



(28) 



(29) 



Substituting the series expansion for f(x) into 7i, we have 



ft oc 

h®* / ^oVC'E a 

00 

= (7oV4)- 1 E 

n— — c 

00 

-(7o 2 /4)- 1 E 



e -(2<n £ + 7 )(i-T) dT 



e inn m t _ e -(7+2ini)t 

7 + « (2£1l + nfi m ) 



n— — 00 



7 + z (2Ql + nQ, m ) 



(30) 



where we have discarded the transient term e (7+2«nt)* < 
The expression for /2 [Eq. I|26|) ] can be found from that 
forli: 



hit) 



(7g/4)- x E 



6„e 



(31) 



7 + i (2f2 £ + nfi m ) 



where the coefficient b n is defined by 

da„ 



-A 



2 ^ 

M 2 



(32) 



In order to find the relative values of a n and b n , it is useful 
to have an approximate expression for them. Assuming 
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that 70 <C Ao, we can replace f(x) with a delta function 
normalized to the same area: 



7T7Q 



S(Dq — cosx) . (33) 



Substituting this expression into Eq. 1)29(1 . we obtain 
7o cos (n arccos D ) 



2A v^~D2 



(34) 



This approximation breaks down for \Dq\ within ~7o/Ao 
of unity. However, as we see below, we are interested 
in integrals of a n over effective detuning, which can be 
well approximated using the expression (f 34|> . We are also 
limited by this approximation to harmonics n <C Ao/70, 
since the factor cos nx is assumed to not vary rapidly 
over the optical resonance. Thus, from Eq. (|32|1 . b„ can 
be approximated by 



70 cos (n arccos Dq) 



4A 
1 



n 2 



(35) 



Thus we see that I2 — I\/2 and the terms of Eq. I|24|) 
proportional to cancel. Substituting Eq. lj3U|) into 
Eq. TMty . we obtain 



4il 2 N ^ 
^21 - ^7— 2^ 



a n e 



37o 



n— — 00 



(36) 



The result (|36|l applies to atoms that are at rest. We 
now modify this result to describe an atomic ensem- 
ble with a Maxwellian velocity distribution leading to 
a Doppler width of the transition. For an atomic 
velocity group with component of velocity v along the 
light propagation direction, the light frequency is shifted 
according to u>(v) = u> (1 — v/c) = lj — kv where k is 
the light-field wave number. Writing the dimensionless 
Dopplcr-shift parameter D v — —kv/Aa, the atomic den- 
sity N for this velocity group becomes 



dN(v) = 



Ap 



Noe -(D v A /r D r dDvj (37) 



where Nq is the total atomic density, and the average de- 
tuning parameter Dq becomes Dq(v) = Dq + D v . Defin- 
ing the velocity-dependent coefficient a n (v) by 

a n (v) dD v = — a n 

70 jV 

_ cos [n arc cos (D + D v )\ c - (DvAn /r D ) 2 dD 



(38) 



2^Jl-(D + D v y 



the velocity-dependent ground-state coherence p2i(v) is 



given by 



p 2 i(v) dD v 



4ft 2 dN(v) 
37o 

4fl 2 N 
3T D 



^ 7 + i (2Q L + nCl m ) 

n— — 00 ' 



E 



a n (v) e 



n— — 00 



7 + i (20, L + n£l m ) 



dD r . 



(39) 



In a situation in which atomic collisions are important, 
such as in a vapor cell with a buffer gas or an antire- 
laxation coating, this result must be further modified to 
take into account collisionally induced velocity mixing. 
For atoms contained in an antirelaxation-coated vapor 
cell, we assume that each velocity group interacts sep- 
arately with the excitation light, but after pumping all 
groups are completely mixed. This model applies when 
light power is low enough so that optical pumping av- 
eraged over the atomic velocity distribution and the cell 
volume is unsaturated. The ground-state coherence of 
each velocity group becomes the velocity-averaged quan- 
tity f>2i(v), given by the normalized velocity average of 
Eq. 

dN(v) f°° . . , 
p 2 i(v)dD v ~—- — / p 2 i(v)dD v 



No 

m 2 dN( v ) 

3Td 



V 



71— — OO 



7 + i (2Cl L + nft m ) ' 

(40) 

where the averaged coefficient a n is given by 
a n (v) dD v 

cos [n arccos (Dp + D v )} p -( DvAo /r f dD 

(41) 



2Ml-(D + D v 



Below, we will need the real and imaginary parts of P21, 
given by 



Rep 2 i(v) dD v 



AO 2 dN(v) 
3Tb 



a n [7 cos nft m t + (2Ql + Ti^m) sinnQ m t] 



n— — oo 



7 2 + (2fi L + n n m y 



lmp2i(v) dD v 



4il 2 dN{v) 
3Tb 



E 

n=—oo 



a n [7 sin nVlmt — (2il^ + nft m ) cosntt m t] 



7 2 + (2n L + n n m y 



(42) 



G. Optical properties of the medium 



We now derive the formula for the optical rotation in 
terms of the polarization of the medium P = Tr pd. The 
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electric field of coherent light of arbitrary polarization 
can be described by [23| 



E : 



i r 

2 
i 



P o e i0 (cos ip cos e - i sin ip sin e) e i{u}t kz) + c.c. 
Eoe^ (sin ip cos e + i cos 99 sin e) ^M-^) + c .c. 



(43) 



where k = lu/c is the vacuum wave number, <fi is the 
overall phase, ip is the polarization angle, and e is the 
ellipticity. 

Substituting Eq. (|43|l into the wave equation 



dz 2 c 2 dt 2 



E 



4tt d 2 . 



(44) 



and neglecting terms involving second-order derivatives 
and products of first-order derivatives (thus assuming 
that changes in <p, e, and <fi and fractional changes in 
Eq are small), gives the rotation, phase shift, absorption, 
and change of ellipticity per unit distance: 



sec 2e [cos ip (Pi sin e + P4 cos e) 

— sin ip (P 2 cos e — P3 sin e)] 



dip 


2iru> 


dz 




# 


2ttlo 


dz 




dEo 

dz 


2ttuj 
= Ism 

c 


de 

dz 


2ttlo 
= 77 sin 



sin y> (P 2 sin e — P 3 cos e)] , 
P4 cos e) 

cos tp (P 2 cos e + P3 sin e)] , 



cos ip (P 2 sin e — P3 cos e)] , 
(45) 



where the components Pi, 2, 3,4 of the polarization are de- 
fined by 



P = i 

2 

1 

+ 2 



(Pi - «P 2 ) e l{ult ~ kz) + c.c. 
(P 3 - iP A ) e'l 1 "" 12 ) + c.c. 



(46) 



For initial values of ip = e = 0, the rotation per unit 
length is given by 



dip 2tt L0P4 

~dl ~ cPn 



(47) 



H. Calculation of the signal 

We now evaluate P = Trpd and substitute into Eq. 
(|47|l to find the optical rotation in terms of the ground- 
state atomic coherence derived above. Taking into ac- 
count that in the nonrotating frame the optical atomic 
coherences oscillate at the light frequency to, we find for 
the polarization components 



P l = V 3 (£,gJg\\dUeJe) Re (p 31 - p 23 ) 



Pi = -\\ 3 Ks-^IMIICe-/e)Im(p31 +/?23) • 

P.-i = = -\/|(4.,^/ii f/ !!v. J . > Til) (A-ij /'2:;i . 
2 



(48) 



^4 = -\/o (^ J fflMll^eJe)Re(/331 + P23) • 



so the optical rotation angle per unit length is given by 



dip _ TTUJ {£gJg\\d\\£ e J e y 

de ~ 3Qc 

7oA 2 



Re (p 3 i + p 32 ) 



(49) 



167rft 



Re (/O31 + P32) 



where A is the transition wavelength. Here we have used 
the fact that for a closed J —> J' transition |21| , 



(50) 



and that u> ~ ojq. 

Substituting in the expressions l|2(J|) . and assuming 
pn ~ p 22 ~ AT/ 3 and S1l <C 70, Eq. l|4^|l can be written 
in terms of the ground-state coherence as 



dp 7o A 2 / Q L (N/3 - Rep 2 i) + (70/2) Impai 2U L A 2 (N/3 - Rep 21 ) \ 

d£ Stt ^ 7o 2 /4 + A 2 (7o 2 /4 + A 2 ) 2 J' [ ' 

or, for the case of complete velocity mixing: 

dtp(y) ,p _ 7qA 2 / nL[djV(«)/3-Rep2i(u)d£> 1) ] + (7o/2)Im J 02i(«)d£>», 2tl L A 2 (v) [dN(v)/3 - Rep 2 i(v) dD v ] \ 
dl " 8n { 7o74 + A» [ 7o 2 /4 + A 2 (z;)] 2 )' 

. (52) 



where the velocity-dependent effective detuning A(v) is given, as before, by 

A(v) = A [D (v) - cos fi m r] . (53) 
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The in-phase and quadrature signals (see Sec. [HJ per 
unit length of the medium, measured for a time T at the 
j-th harmonic of the modulation frequency, are given by 
the time averages 



df) n (v) dD v 
■ dD„ 



df 



d<p(v) 



T J dl 



cos (J f2 m £) dt, 



dtp^iv) dD v 
■ dD„ 



df 



T 



T dip(v) . , . _ „ , 
-ZjU sin j fl m t) dt. 
df 



(54) 



We substitute the formulas for the real and imaginary 
parts of the ground-state coherence [Eq. (|42|) ] into the 
formula for the optical rotation [Eq. i|52|l ]. and the re- 
sulting expression into Eq. I)54[l . After evaluating the 
time integrals (see Appendix |A"|) . we find that the signals 
due to each velocity group are given by 



df 



dip° ut (v) 
df 



dD v 



dD v ~ - 



A 2 » 2 
67rr D 

oo 

E 

n— — oo 

A 2 n 2 



7 2 + (2fi L +nr! m )^ 



67rrc 



di\r(u) 

OO _ , i 

\ -« 7 &n l^n+j — a n-j ) 

„±^,7 2 + (20l + nfi m ) 2 



(55) 



Using the definitions of dN(v) and a„(v) [Eqs. (|37I38[) ] 
we can rewrite Eq. I|55|l as 



dD„ 



E 



7 (2f2 L + nVL m ) a n [a n +j(v) + a n -j(v)] 



1 2 + (2n L + n n m y 



dD v , 



dip° ut (v) 
d£ 



dD v 



V E 



n— — oo 



7 2 a n [a n +j(v) - a n -j(v)] 
7 2 + {2Q L + nfl m ) 2 



dD v , 



(56) 



where the signal amplitude factor rj is defined by 



V = 



1 » 2 7o 

6tt r 2 7 



A 2 iV . 



(57) 



The total signal, given by the integral over all velocity 



groups, is found by replacing a n (v) with a n 



df 



dip° ut 



total 



dipin(v) 
df 



dD„ 



df 



OO 

^E 

n— — oo 

/■OO 



7 {20, l + ^m) a n (a n+ j + a n -j) 



7 2 + (20 ^ + nO ri 



total 



d(p ou t{v) 

dt 



dD v 



7 2 a n (a n+ j ~ a n -j) 
n ^ oc l 2 + (2n L +nQ m ) 2 



E 



(58) 



Each term of the sums corresponds to a resonance at 
f2k/fl m = — n/2 fFies. 13*14^1 . Near each resonance the in- 
phase signal is dispersive in shape, whereas the quadra- 
ture signal is a Lorentzian. When plotted as a function 
of the Larmor frequency normalized to the modulation 
frequency, f2z,/f2 m , the widths of the resonances are de- 
termined by the normalized ground-state relaxation rate 
7/fi m . The relative amplitudes of the resonances are de- 
termined by the ratio of the modulation depth to the 
Doppler width, Ao/Id, and the normalized average de- 
tuning A;/A . 



IV. CONCLUSION 

We have presented a theory of nonlinear magneto- 
optical rotation with low-power frequency-modulated 
light for a low-angular-momentum system. The 
magnetic-field dependence predicted by this theory is 
in qualitative agreement with experimental data taken 
on the Rb Dl line. Directions for future work in- 
clude a more complete theory describing higher-angular- 
momentum systems, including systems with hyperfine 
structure, and higher light powers. A possible compli- 
cation to the FM NMOR technique in systems with hy- 
perfine structure is the nonlinear Zeeman effect present 
at higher magnetic fields, so a theoretical description of 
this effect would also be helpful. FM NMOR has been 
shown to be a useful technique for the selective study of 
higher-order polarization moments, which give rise to dis- 
tinct resonances at different values of the magnetic field 
than the quadrupole resonances studied here ja] (see also 
Ref. Higher-order moments are of interest in part 

because signals due to the highest-order moments possi- 
ble in a given system would be free of the complications 
due to the nonlinear Zeeman effect. To describe these 
moments, a calculation along the same lines as the one 
presented here but carried out to higher order and in- 
volving more atomic sublevels would be necessary. 
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APPENDIX A: EVALUATION OF THE TIME 
INTEGRALS 

In evaluating Eq. I|54(l . several time integrals appear: 

rfi m T 



h 



1 



Q T 



1 



fi T 



h = 



l 



Q T 



cos JX 


7o 2 /4- 


^A 2 [D (v) - 


i 2 
cos x] 




sin jx 




7 74- 


^A 2 [D (v)- 


cos x] 




cos jx sin nx 




7 2 /4 - 


l~ Aq [D (v)- 


cos x] 




sin jxsinnx 




7 2 /4 ^ 


l~ Aq [D (v)~ 


cos x] 




cos jx cos nx 




7o 2 /4- 


l~ A 2 , [A)(V>- 


l 2 
cos x] 




sin jx cos 





dx, 



dx, 



dx, 



dx, 



o 7o 2 /4 + A 2 , [£> («)- cos x]' 



dx, 



as well as the related integrals 



h = 



Aq [£>o( w ) — cosx] 2 cosjx 
{ 7 o 2 /4 + A§ [^ ( v )-cosx] 2 }^ 



■ dx 



= -A 



°SA 2: 



Aq [Do(v) — cosx] 2 sin jx 
{ 7o 2 /4-t-A 2 [D (v)-co S x} 2 Y 



■ dx 



hi 



Q T 



7H-L JQ 



Aq [Dq(v) — cosx] cosjxsinnx 
{ 7o 2 /4 + A 2 [D (v)-co S x} 2 } 2 



dx 



2 9/ 5 

°dA 2 ' 



/is 



1 



O T 



Sl m T A 2 



Aq [Do(v) — cosx] sinjxsinnx 
{ 7o 2 /4 + A2 [A>( v )-cosx] 2 } 2 



dx 



-A 



g/e 

dA 2 ' 



1 



'13 



n m T a 2 



"m/ 1 JO 



_ A 2 



Aq — cos z] cos jx cos n. 

{ 7o 74-t-A 2 [^o( v )-cosx] 2 }' 



■ dx 



9A 2 ' 



'14 



Jo 



A 2 



Aq [_Dq(u) — cosx] sinjxcosnx 
{ 7 2 /4 + A 2 [iJ ( V )-cosx] 2 } 2 



dx 



'9A 2 ' 



If T is many modulation periods, the above integrals can 
be approximated by averages over one period. Thus, we 
can change the limits of the integrals to (— tt,tt), and set 
the normalizing factor (O m T) to (2tt) . Using the 
trigonometric substitutions 



cos jx cos nx 
cos jx sin nx 
sin jx sin nx 
sin jx cos nx 



{cos [(n — j) x] + cos [(n + j) x]} , 
{sin [(n - j) x] + sin [(ra + j) x]} , 
{cos [(n — j) x] — cos [(n + j) x]} , 
^ {sin [(n - j) x] + sin [(n + j) x]} , 



(Al) 
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we can rewrite the above integrals in terms of the a n and 
b„ coefficients. 



h 




-l 

CLi, 




L 

-*-4 


~ 






h 


~0, 








— W0/ 4 J 


-l (a n -j 


— O-n+'j ) 

2 


h 


- (7o 2 /4) 


-i (a n -j 


+ a n+j ) 
2 


h 








h 


(7n 2 /4) 

V 101 J 


bi. 




Jin 
± 10 


~ 






ill 


~0, 






h2 


- (7oV4) 


-1 (b n -j 


- b n +j) 

2 


h 3 


- (7o 2 /4) 


-1 {bn-j 


2 


hi 


~ 0. 







(A2) 



As in the evaluation of Eq. I|24|l . use of the approximate 
expression b n — a n /2 [Eq. ffi5jl ] results in the cancella- 



tion of some terms proportional to Ql, producing the 
relatively simple form of Eq. H55|) . 
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